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Fundamental measure theory for mixtures of parallel hard cubes. 
Phase behaviour of the one-component fluid and of the binary mixture 
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A previously developed fundamental measure functional [J. Chem. Phys. 107, 6379 (1997)] is 
used to study the phase behaviour of a system of parallel hard cubes. The single-component fluid 
exhibits a continuous transition to a solid with an anomalously large density of vacancies. The 
binary mixture has a demixing transition for edge- length ratios below 0.1. Freezing in this mixture 
reveals that at least the phase rich in large cubes always lies in the region where the uniform fluid 
is unstable, hence suggesting a fluid-solid phase separation. A method is developed to study very 
asymmetric binary mixtures by taking the limit of zero size ratio (scaling density and fugacity of 
the solvent as appropriate) in the semi-grand ensemble where the chemical potential of the solvent 
is fixed. With this procedure the mixture is exactly mapped onto a one-component fluid of parallel 
adhesive hard cubes. At any density and solvent fugacity the large cubes are shown to collapse into 
a close-packed solid. Nevertheless the phase diagram contains a large metastability region with fluid 
and solid phases. Upon introduction of a slight polydispersity in the large cubes the system shows 
the typical phase diagram of a fluid with an isostructural solid-solid transition (with the exception 
of a continuous freezing) . Consequences about the phase behaviour of binary mixtures of hard core 
particles are then drawn. 



PACS: 61.20.Gy, 64.75.+g, 82.70.Dd 



I. INTRODUCTION 

This paper is the sequel of a previous one§ (henceforth 
referred to as I) in which the so-called fundamental mea- 
sure theory (FMT) was applied to build a density func- 
tional for the multicomponent system of parallel hard 
cubes (PHC). In I we explained all fundamentals of the 
theory and gave a full account of the technical details 
involved in the derivation of the functional. We also dis- 
cussed the pros and cons of the theory, as compared with 
more standard density functional theories (DFTs), and 
suggested some possible extensions. 

In this paper we apply the formalism developed in I to 
study the phase behaviour of the PHC fluid, with spe- 
cial emphasis in its relevance for the understanding of 
the phase behaviour of mixtures. The PHC fluid is a 
rather academic model and it possesses a bunch of 'pe- 
culiarities' which are rather odd for a fluid model, e.g. 
the uniform fluid is anisotropic at small scales (the cubes 
are kept parallel to each other), freezing occurs at very 
low packing fractions (around 0.3-0.4) and it is a contin- 
uous, instead of first order, transition (a consequence of 
the lack of isotropy of the fluid phase), and the depletion 
in the binary mixture is very strong compared to hard 
spheres (HS). But in spite of these peculiarities — which 
certainly make of this model a caricature of a fluid, the 
physics one can learn from its phase behaviour can be 
easily extended to more reasonable fluids, and it has the 
important added value of being a much simpler model to 
carry out analytical calculations (even more: FMT seems 
to be somehow optimal for this modelEra). 

There have been a few pgewious studies in the litera- 
ture about the fluid of PHCHeI scattered in the last forty 



years, but the relevance of this model has only recently 
become apparent when it has been proposed as a model p£ 
a fluid able to demix by a purely entropic mechanism Jj€2I 

Entropic demixing has been a long standing question 
which only recently begins to be understood. It is well 
known that different attractions between two types of 
particles in a fluid can produce segregation into two 
phases, each rich in one type of particles .lid The ques- 
tion remains whether hard particles, for which only an 
entropic balance can drive a phase transition, ever demix 
and how. It is clear that nonadditive mixtures (mix- 
tures in which particles of diffcrentj-typfi interact as if 
they had a larger volume) do demix,E3~t3 but they have 
the segregation mechanism introduced at the interac- 
tion potential. So the nontrivial question concerns ad- 
ditive hard-particle mixtures. The question is tricky be- 
cause the simplest model of this type=pHS — was solved 
in Percus-pYjevick (PY) approximatiorJIS and shown never 
to demix.EZI It had to wait almost thirty years until the 
PY result was questioned. By solving numerically the 
Ornstein-Zernike equation with closure relations more ac- 
curate than the PY one a spinodal instability was shown 
to occur in a binary mixture of HSEjEj for a diameter 
ratio below 0.1-0.25 (depending of the authors). 

It was then believed that a sufficiently asymmetric HS 
binary mixture undergoes fluid-fluid demixing in a cer- 
tain region of the phase diagram. But successive exper- 
iments performed in suspensions of polystyrene or silica 
spheres (which to a large accuracy can be considered HS) 
showed cumulative evidence that demixing is coupled to 
freezing, and so one of the phases shows. up, as a crystal — 
sometimes a glass — of large particles Oj£jSJlecent theo- 
retical calculations confirm this scheme. c§0 
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Direct computer simulations have run into troubles 
when dealing with this problem due to the extremely low 
probability of moving a larger particle in a sea of small 
ones without overlapping any of them. Simulations are 
still possible for not top.dissimilar diameters, but they 
do not show demixing .E3E3 With the help of specially 
designed cluster moves the first evidence of demixing has 
recently been found in simulations for diameter ratios 
smaller than 0.05.E3 Unfortunately these cluster moves 
do not help above the percolation threshold, and this sets 
a relatively low upper bound for the total packing frac- 
tion of the fluid. Besides, this method does not allow to 
identify the coexisting phases, so it does not distinguish 
fluid-fluid or fluid-solid demixing. 

In order to elucidate the nature of demixing for very 
asymmetric mixtures the attention has shifted to under- 
stand the depletion interaction. An effective pair poten- 
tial between theJ|afae spheres can be obtained by differ- 
ent procedures.c2l~E3 Its shape reveals a very deep and 
narrow (one small-sphere diameter) well followed by a 
couple of small oscillations extending two or three small- 
sphere diameters. The depth and the amplitude of the 
oscillations depend on the small spheres packing fraction. 
In view of the phase behaviou r . of . spheres with narrow 
and deep attractive potentialscaH fluid-fluid demixing 
is ruled out for sufficient asymmetry, and instead fluid- 
solid demixing or even expanded-dense solid demixing 
(i.e. demixing once the large component has crystallised) 
should appear. Simulations of a.-syptem of HS supple- 
mented by the depletion potentiaErEj confirm this con- 
jecture if the diameter ratio is below 0.1. Expanded- 
dense solid demixing isindecd shown in some of this sim- 
ulations for ratios 0.1l3 or 0.05a This phase behaviour 
has been recently corroborated by direct simulations of 
the binary mixtureJ23 which have been shown to be pos- 
sible in the relevant region of the phase diagram thanks 
to the small amount of small spheres present in those 
statepoints. 

In the limit of zero size of the small spheres at constant 
packing fraction th e dep letion potential becomes Baxter's 
adhesive potential.EZrEa By applying a similar limit to the 
FMT functional of a binary mixture of PHQm we have 
recently obtained a functional for parallel adhesive hard 
cubes (PAHC).l3 By avoiding the singularities of this po- 
tential (we will treat this point in detail later on) we show 
that the phase behaviour of this fluid is consistent with 
the above described picture for the mixture of HS, with 
the difference that expanded-dense solid demixing is the 
most common scenario because PAHC freezing occurs at 
rather low packing fractions. 



II. THE ONE-COMPONENT FLUID 

There are relatively few results in the literature con- 
cerning the fluid of PHC. As concerns the uniform fluid 
the virial coefficients, both in 2D and 3D, have been ex- 



actly obtained through diagrammatic expansions up to 
the seventh.Ero They have also been obtaincjd for differ- 
ent approximate integral equation theories.B The main 
consequence one draws from this information is that the 
virial expansion is very poorly convergent. For the 3D 
fluid, the seventh order virial expansion exhibits a maxi- 
mum in the pressure at 77 f=a 0.6, .and then goes down very 
quickly to reach negative values!] Its behaviour at moder- 
ate densities strongly deviates from that of a sixth order 
virial expansion — in contrast with what happens for HS. 
This is the fingerprint of a nearby divergence. In fact, the 
simulations of the fluid of PHC (3D)EH3 [ ahow a continu- 
ous freezing into a simple cubic lattice at£l rj = 0.48 ±0.02 
(in contrast with the first order nature of the HS freez- 
ing). This result has been proven to be exact in infi- 
nite dimensions a The reason for a continuous freezing is 
the lack of rotational symmetry of this system even in 
the disordered phasea (thus freezing does not brakes this 
symmetry). Allowing the cubes to rotate^restores the 
first order nature of the freezing transition.!] 

FMT's equation of state of the uniform PHC fluid is 
simply that of the scaled particle theory (SPT), i.e., from 
Eqs. (60) and (61) of I, 

PP/ P =J--^L, (D = 3), (2) 

with /3 the inverse absolute temperature in Boltzmann 
constant units, P the pressure, p the number density, 
and 77 the packing fraction, r\ = pa D , a being the cube 
edge-length. In Fig. |l| the SPT equation of state @ is 
compared with the simulation data. 

On the other hand freezing can be studied as usual in 
DFT by parametrising the local densiti£-,with a sum of 
gaussians centered at the lattice points.c3 Since the lat- 
tice is simple cubic, this can be easily achieved by setting 

00 

p ^ = n n ^ ~ n " rf ) ' ( 3a ) 

u=x,y,z n l/ —~oo 
/ ry \ 1/2 2 

g{u) = (-) , (3b) 

with d = r\~ x l^(j the lattice spacing, and a. a variational 
parameter determining the localisation of particles at the 
lattice sites. 

In I we expressed the free-energy functional of a multi- 
component system with density profiles Piijc) (i labeling 
the species) as 

F[{ Pi (r)}\ = F id [{ Pt (r)}} + F™{{p t (r)}} , (4) 

/3F id [{p 4 (r)}] = W dv Pl (r) [lnV iPi (r) - 1] , (5) 

i J 

pF e *[{p i (r)}}= J dv<t>({n a (r)}), (6) 
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with Vi the thermal volume of species i, and n a = 
J2i Pi * w i a se t of four weighted densities (stars de- 
notes convolution) defined by the weights 

4 0) =6?5V5*, (7a) 

wf'^W.W.W), (7b) 

wf = (W?,^d«flf,ef^f), (7c) 
= efept , (7d) 

with 0f = 9(cTi/2 - M), and Sf = (1/2)5(^/2 - |«|). 
As in I we will also need the two scalar densities = 
Ui t x + rii,y + rii tZ , rii yVl v = x, y, z, being the vector com- 
ponents of n i; i — 1,2. The function $(r) is given in 
terms of the n a 's as (see I) 



$(r) = -n ln(l-n 3 ) + 



1 - n 3 (1 - n 3 y 



with U2,!/, u = x,y, z, the components of ri2. 
For convenience let us introduce the functions 

e(u) = -erf (^/auj , 



(8) 



(9) 



p(u) = ^2 [e(u-nd + a/2) - e(u-nd- a/2)], (10) 

n— — OO 
OO 

q(u) = [g(u~nd + a/2) + g(u-nd- a/2)} (11) 



n— — oo 



[notice that g(u) — e'(u)]. The latter two are periodic 
with period d. In terms of these functions 

n 3 = p(x)p(y)p(z) , (12a) 
n 2 = ^{q(x)p(y)p(z),p(x)q(y)p(z),p(x)p(y)q(z)) , (12b) 

n i = k{p(x)q(y)q(z),q(x)p(y)q(z),q(x)q(y)p(z)) , (12c) 



"0 = ^q(x)q(y)q(z) , 

and accordingly <&(r) = no(r)^(r), where 
^ = - ln(l - n 3 ) 
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(12d) 



(13) 



1 - n 3 (1 - n 3 ) 2 
This simplification occurs because of the factorisation of 



the density ( |3a| ) — and hence of (12). 

The free energy per particle, \& = (JF/N, of the solid 
phase can be obtained as the integral over a unit cell of 
the lattice of the free energy density. Hence the ideal 
part turns out to be 



\I/ id = 



rd/2 


fd/2 


fd/2 


/ dx 


1 dy 


\ dzp{r)\\nVp{v)-l] 


l-d/2 J 


-d/2 J 


-d/2 



ln( V/a 3 
+3 



/OO / °° \ 

dx g(x) In I crg(a; — nd) J 

-°° \n=-oc / 



(14) 



The integrand is written in a suitable way to use Gauss- 
Hermite numerical quadratures .cil As for the excess con- 
tribution, similar manipulations lead to 



^ rd/2 f-d/2 ^d/2 

* ox = -/ dxq(x) dyq(y) dzq(z)ip(x,y,z) 

8 J — d/2 J -d/2 J -d/2 

/>oo />oo y* oo 

dxg{x) I dyg{y) \ dzg(z) 



x^2ip(x± a/2, y ± a/2, z ± ct/2) , 
{±} 



(15) 



where the last summation runs over all combinations of 
signs. Again the latter expression is suitable for Gauss- 
Hermite integration, which is crucial this time because 
dl5| ) involves a three-dimensional integration. 

We can now minimise with respect to a to determine 
the equilibrium profile. This yields a continuous freezing 
at rj = 0.348. As the transition is continuous we can make 
a more accurate determination of the transition density 
via a standard bifurcation analysis. This is equivalent to 
finding the density at which the structure factor diverges 
for some wavevector k c . The structure factor is expressed 
in terms of the DCF as 



5(k) 



1 - pc(k) ' 



(16) 



and the Fourier transform of the DCF, c(k), is obtained 
through Eq. (56) of I for a one-component fluid. In order 
to simplify the final expression we exploit the symmetry 
of the crystal by choosing k c = (fc c , 0, 0) (the result would 
be the same if we chose k c along the Y or Z axes) . Thus 
the condition to determine the critical point is 



4-3?/ c + r;2 . 
1 + 2? ?c f1 — Jo(ka) 



n, 



(i - Vc) 3 
a 9 - 47fc + j]l 



Jo(W2) 2 >0, 



(17) 



the equality holding only for k — k c ; jo(x) = smx/x is 
the zeroth order spherical Bessel function. The solution 
to this equation is r\ c = 0.3143 . . . and k c a — 4.8276 .... 

It is noticeable the discrepancy between the value of 
•q c obtained from the divergence of the structure factor 
and that obtained using the profile (|). The reason for 
this discrepancy can be inferred if we obtain from k c 
the lattice spacing, d c /a — 2ir/k c — 1.3015...; thus 
"Hc{d c /a) 3 = 0.6929 . . . , what means that the resulting 
crystal has a large fraction of vacancies (around 31% of 
the lattice sites!). This is a strong effect that can be ac- 
counted for by simply multiplying the r.h.s. of ( |3aj ) by an 
average occupancy ratio § and minimising with respect to 
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this new variational parameter. In the calculation pro- 
cess this simply amounts to (i) add a term lni? to the 
ideal free energy (| l4T) : and (ii) replace ns by $713 in the 
definition of ip [Eq. Q13)] . Notice that d is now given by 
d = a^/ri) 1 / 3 . As a result we obtain the correct value of 
r) c and an occupancy ratia-pf ~ 0.694, consistent with 
the value obtained above.c3 



the matrix 




FIG. 1. Equation of state of the PHC fluid. Solid lines 
correspond to the pressure of the stable phase (fluid or solid) 
at the given packing fraction, 77. Dotted line is the unstable 
fluid branch beyond freezing. Dashed line is the solid branch 
computed without accounting for vacancies in the lattice (see 
text). Full circles are the simulations of Ref. ^|and empty cir- 
cles those of Ref. [?] (actually we have averaged the two data 
sets reported for the largest system size). 

The solid equation of state — with and without 
vacancies — is plotted in Fig. |l| and compared with the 
simulations. It is obtained as /3Per 3 = r\ 2 d^ldr\. We 
can see that the overall agreement is good, although the 
freezing point is shifted down with respect to the simu- 
lations because the SPT equation of state overestimates 
the pressure of the fluid phase. 



III. STABILITY OF THE BINARY FLUID 
MIXTURE 

Two— are the requirements for a fluid mixture to be 
stablerJ (i) the positiveness of the specific heat at con- 
stant volume (cy), and (ii) the positive definiteness of 



M, 



d 2 f 
dpidp-j 



(18) 



where / = F/V, F being the Helmholtz free energy, and 
Pi the number density of species i. Condition (i) is triv- 
ially fulfilled, because any hard core model is athermal, 
what means that the dependence of the free energy on 
temperature is that of an ideal gas; hence the positiveness 
of cy . Condition (ii) is a consequence of the equilibrium 
state being a minimum of the free energy. 

According to M) the matrix M will be given by 



X>. w 



dn-ydria ' 



and relating the second term of the r.h.s. of this equa- 
tion to the DCF of the fluid [Eq. (7) of I] it can simply 
be written as 



M, 



— Sij - Cij(0) . 
Pi 



(19) 



Now, for a binary mixture, M is a 2 x 2 matrix with 
all its elements positive; thus the mixture will be stable 
provided |M| > 0. The solution to the equation |M| = 0, 
if it exist at all, will represent a spinodal curve. Such 
a condition can be understood in terms of the structure 
factor matrix of the mixture, given by 



pstk^p-i-cfk)]- 1 , 



(20) 



where Q 3 = Cij and = /0j<%. Then, after (|Ts|), 
|M| = is the condition for the structure factor to di- 
verge at zero wavevector (the uniform fluid). 

In order to work out the expression of |M| let us in- 
troduce the following notation: r\i = af pi, the packing 
fraction of species i; r\ = rji + 772 = £,d, the total packing 
fraction of the fluid; r = CT1/CT2, the large-to-small edge 
ratio (r > 1); and x = rji/rj, the relative packing fraction 
of the large component. With these definitions as well as 
the short-hand M. = P1P2IMI we obtain, after a tedious 
but straightforward calculation, 



M 



1 + 7? 

(1 - 77)3 



(21) 



for D — 2, an 



M 



i 4 1 3(r-l) 2 

1 + - + - — -x(l - x) 

77 i] z r 



(22) 



for D = 3. 

From ( ^l| ) it follows that the 2D mixture is stable 
whichever the values of rj, r and x. Accordingly, parallel 
hard squares never demix into two fluids with different 
composition. Equation (|22|), however, tells us that the 
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mixture of PHC will be stable provided the expression in 
square brackets is positive, i.e. 



where y = r)/(l — rj), and /3A/ij = fim - ln(Aj/cr; 



1 



4 
V 



1 

if 2 



> 



3(r-l) : 



x(l — x) . 



(23) 



6 > 



- lOr + 1< . 



Since the minimum of the function 1 + 4/rj + I ft] 2 , for 

< t) < 1, is 6 (reached when rj = 1) and the maximum 
of x(l — x), for < x < 1, is 1/4 (reached when x = 1/2), 
( ^3| ) will hold for any n and x whenever 

3(r - l) 2 1 
r 4 

which means for any l/r c < r < r c , where 

r c = 5 + VM « 9.98 (24) 

(or, equivalently, as we have defined r > 1, for any 

1 < r < r c ). For r > r c there will exist values of rj 
and x for which (23) does not hold, and thus the mixture 
demixes. From (23) it is very simple to find that those 
values correspond to the region above the curve 



(r-iy 



x(l — x) 



(25) 



which therefore defines the spinodal. Figure ^ shows this 
curve for a few values of r. It is interesting to notice the 
symmetry of the spinodal with respect to x = 1/2. This 
means that for a given packing fraction, rj, the stability 
of the mixture depends on the fraction of occupied vol- 
ume of any of the particles, regardless their type, large 
or small. 

The existence of a spinodal instability of this type 
means that if the system is kept at constant pressure, 
there is a certain region in the density-composition phase 
diagram in which the fluid is stable in two coexisting 
phases, one rich in small particles and the other one rich 
in large particles. In order to determine the values of rj 
and x of the two coexisting phases we must solve, at a 
given pressure, the equilibrium equations. If we denote 
rj s , x s (respectively rji, x{) the values of the small-particle 
(respectively large-particle) rich phase, those equations 
can be written 

P(jj s ,x s ) = P(r)i,xi) = p, 
(j,x(r) a ,x a ) = Hi(r/i,Xi) , 

p being the externally fixed pressure. This equations 
express equality of the pressure and chemical potentials 
of both kinds of particles in each of the two coexisting 
phases. For the present binary mixture, jXi, (12 and P 
are given by 

f3Pol = y(r 3 - (r 3 - l)x) + 2y 3 (r - (r - 1» 3 

+3y 2 (r 2 - (r 2 -l)x){r - (r -l)x) , (26) 

(3Afi! = ln(xy) + 3y(r 2 + r - (r 2 + r - 2)x) 

+3y 2 (r-(r~l)x) 2 +{3Pa 3 1 , (27) 



/3A/x 2 = ln[(l - x)y] + 3y(2 - (2 - r" 1 - r~ 2 )x) 



S=r 0.5 




+3^(1 - (1 - r- L )xY + pPcrfr 



3--3 



(28) 



1.0 



FIG. 2. Phase diagram of the demixing transition of par- 
allel hard cubes; n is the fraction of volume occupied by all 
cubes, whereas x = Vi/v is the fraction of volume occupied 
by the large cubes. Short-dashed lines represent the spinodals 
for different values of the edge-to-edge ratio; solid lines are the 
corresponding coexistence lines (the actual transition lines); 
the dotted line is the line of critical points of the demixing 
transition for all values of r > r c ; finally, the long-dashed 
segments joining black dots are two examples of coexisting 
states. 



To determine the critical point of this transition we 
use the fact that this point is the only one for which the 
spinodal and the coexistence line coincide, and it cor- 
responds to the minimum value of the pressure on the 
spinodal. Hence Eqs. ( p5| ) and (|2^) allow us to deter- 
mine this point for every value of the edge-ratio r. 

The phase diagram (Fig. ||) shows a few features that 
are worth noticing. First of all, it is interesting to see 
that the critical line collides, when r —> 00, with the 
x = 1 edge of the phase diagram at a nonzero packing 
fraction. What this suggests is that in this limit the pack- 
ing fraction of the small component goes to zero but it 
still remains a residual depletion between the large cubes. 
This depletion forces the one-component effective fluid of 
large cubes to collapse beyond a certain packing fraction. 
We will explore this matter in full detail in Sec. 0. 

Another interesting feature of the phase diagram is 
the remaining impurity of the two separated phases even 
when the system undergoes an infinite pressure. This 
reflects in the fact that the coexistence lines end up at 
values of x other than 1 or (pure components) when 
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r/ = 1 (actually, this effect is noticeable only for the values 
of the large-cube rich phase, although it is also present 
in the other phase). The prominent asymmetry of the 
coexistence line is another striking feature, but easy to 
understand: it arises from the enormous volume differ- 
ence between large and small cubes necessary to produce 
demixing (notice that demixing begins for r rs 10, and 
this means that large cubes occupy a volume 1000 times 
larger than the small ones). This forces the large-cube 
impurities in the small-cube rich phase to be in an ex- 
tremely low concentration. The wide metastability re- 
gion in the small-cube rich side means that for those 
compositions the mixture is less sensitive to variations 
in composition. 



IV. FREEZING OF THE BINARY MIXTURE 



In order to check to which extent the demixing scenario 
found in the previous section holds we have to determine 
whether the fluids are stable against spatial modulations 
at the coexisting compositions. Spatial inhomogeneitics 
cause a divergence of the structure factor matrix ( pp| ) at 
a certain nonzero wavevector. Thus for a given composi- 
tion, x, the spatial instability is found as the lowest total 
packing fraction at which the determinant 




M(k) 



-C(k) 



(29) 



0.0 0.2 0.4 0.6 0.8 1.0 
X 

FIG. 3. Phase diagram of the demixing transition of par- 
allel hard cubes; r\ is the fraction of volume occupied by all 
cubes, whereas x = rji/rj is the fraction of volume occupied by 
the large cubes. Dashed lines represent the coexistence lines 
of the demixing transition for different values of r = o"i/<72; 
solid lines are the freezing spinodals for the same values of 
r; finally, the dotted line is the line of critical points for all 
r > r c . 



V. INFINITELY ASYMMETRIC BINARY 
MIXTURE: PARALLEL ADHESIVE HARD 
CUBES 



vanishes for at least one vector k. We can use the expres- 
sion for the DCF found in I [Eqs. (53), (54), and (57) of 
I] and simplify the problem by simply looking for insta- 
bilities along the three coordinate axes. By symmetry, 
this amounts to take k = (k, 0, 0). 

The value of r\ at which .M(k) = 0, as a function of x, 
is plotted in Fig. ||| for different values of r; these lines of 
instability are compared with the coexistence lines of the 
demixing transition (Fig. |^). It is clear from the figure 
that the critical points of the latter are always in the un- 
stable region; therefore, for any pair of coexisting fluids, 
at least one (the large-cube rich one) is always unstable 
against spatial inhomogeneities. In other words, of the 
two phases in which the system phase separates, the one 
rich in large cubes must always be a solid. Notice, on the 
other hand, that the other phase is also unstable for size 
ratios smaller than r ss 30. One is then tempted to con- 
clude that fluid-fluid demixing is preempted by freezing 
in this system. However coexistence between a large- 
cube rich solid phase and a small-cube rich fluid phase 
may change drastically the compositions of the coexisting 
phases and thus make a fluid-solid demixing more stable 
than just a freezing of the whole system. The only conclu- 
sion we can draw from Fig. P is that the fluid-fluid demix- 
ing transition found in Sec. 



A. The binary mixture as an effective 
one-component fluid 

In order to study the phase behaviour of a very asym- 
metric binary mixture let us first consider the effect on 
the interaction of the big particles induced by the small 
ones (depletion). To this purpose let us use a semi-grand 
ensemble in which the small particles (solvent) are kept 
at constant chemical potential. This is the usual exper- 
imental setup in colloidal suspensions. In this situation 
the element of the structure-factor matrix ( |20| ) corre- 
sponding to the correlations between large particles can 
be considered as the structure factor of an equivalent one- 
component fluid made of the large particles interacting 
via the effective potential induced by the solvent. This 
turns out to be a very useful viewpoint. Let us see how 
this come about. 

The appropriate thermodynamic potential for the 
semi-grand ensemble is obtained through a Legendre 
transformation of the Hclmholtz free energy, namely 



III is always metastable. 



T(/i 2 , [pi]) = F[pi,p 2 ] - H2 J drp 2 {v), 

SF , , . SF CX 

= 7-73 = f 3 In V 2 p 2 (r) + 



<Wr) 



(30) 
(31) 



G 



where Eq. ( |3l| ) provides the equilibrium density of the 
solvent for a given chemical potential p 2 and a solute 
density profile pi(rV thus allowing us to eliminate p 2 ( r ) 
from the r.h.s. of (30). The thermodynamic potential T 
can also be looked at as the Helmholtz free energy func- 
tional of an effective one-component fluid, for which p 2 
is just an external parameter tuning the interaction be- 
tween its particles. Accordingly we can separate out the 
ideal and excess parts, 



T(/z 2 , \p x ]) = /T 1 J dr Pl (r) [In V lPl (r) - 1] 

+ T ex (/i 2 ,[ Pl ]), (32) 

where, upon comparison with @), 

T ex (/x 2 , [Pi]) = 0- 1 J dv P2 {v) [lnV 2 p 2 (r) - 1] 

+ F c *[ Pll p 2 ]-p 2 [dvp 2 (r), (33) 



with p 2 (r) determined by Eq. (pl|). 

Now, the DCF of the effective fluid will be 



c e ff(r, r') = -j3- 



5 2 T C 



(34) 



5 Pl (r)5p 1 (r>) 

The first functional derivative of T cx can be written 
ds < j3 lnV 2 p 2 (s) + 



S P2 (s) M2 j S Pl (r) 
(35) 



SF C * 5F C 



8px(r) 5px(r) 
A new derivative yields 



c off (r,r') = c n (r,r') + ^ ds A(r, s)c 2 i(s, r') , 



(36) 



where we have used (|34|), the usual DCF matrix defini- 
tion Cij (r,s) = —[35 2 F cx /5pi(r)5pj(s), and the shorthand 
A(r,s) = Sp 2 (s)/Sp 1 (r). 

The functional A(r, s) can be readily obtained by de- 
riving (plh with respect to p\ (r) , which leads to 



1 



P2(S 

Then 



■A(r,s) - ci 2 (r,s) - / dt A(r,t)c 22 (t,s) = 0. 

(37) 



A(r,s) = J dtc 12 (r,t)B(t,s), 



(38) 



where -B is the solution to 
1 



dt 



P2(r) 



<5(r-t)-c 22 (r,tUB(t,s) = 5(r-s) 



(39) 



Substitution of (|3q) into (36) finally leads to 



c cff (r, r') = cn(r,r') + // dtds c 12 (r, t)B(t, s)c 2 i (s, r') . 

(4 ) 

This expressions have full generality: it is valid for 
any binary mixture of any kind and in any particular 
phase. But we can get a bit further if we particularise 
for the uniform fluid. In this case the density profiles 
are constant and all the above functions depend on the 
difference of their arguments, because of translational in- 
variance. Then, a Fourier transform of Eq. ( |39"| ) permits 
to obtain explicitly B as 



S(k) 



Pi 



1 - p 2 c 22 (k) 



(41) 



Again Fourier transforming ([ 
yields 

c off (k) = cu(k) + 



3) and using ( |4l"| ) finally 
p 2 ci 2 (k) 2 



1 - /3 2 c 22 (k) 



(42) 



In order to understand the meaning of Eq. ( |42| ) let us 
compute the structure factor of the effective fluid: 



PiSeff(k) 



1 



Pi 1 



c 22 (k) 



Pi 1 



Coff(k) 



P- 1 - C(k) 



(43) 



the resulting expression is but the (1, l)-element of the 
structure- factor matrix pS(k), according to its definition 
©. This was, by the way, the starting point from which 
the effective fluid was defined in Ref. |3^. 

Further insight can be gained if we expand the second 
term of ( f42| ) in powers of p 2 : 

00 

c cff (k) = c u (k) + P2 P 2 l ci 2 (k)[c 22 (k)]"c 21 (k) , 

n=0 (44) 

where it can be explicitly seen that while the first term in 
( |42| ) represent the direct correlation between two solute 
particles via the direct potential between them, the sec- 
ond term accounts for the indirect contributions to this 
correlation due to interaction with one, two, three, etc, 
intermediate solvent particles. Is this effect that accounts 
for depletion in the binary mixture of hard particles. 



B. Depletion in the binary mixture of parallel hard 
cubes 



Let us first compare the effective attraction between 
large particles induced by the small ones (depletion). A 
simple way to achieve this is by computing the work we 
have to make against the system in order to separate two 
big particles further than the diameter of a small one. 
This work will be simply PAV, with P the pressure of 
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the fluid and AV the free volume lost by the small parti- 
cles due to the disappearance of the overlap between the 
excluded regions of the large particles (shaded in Fig. |^). 

In the case of HS AV = (3v /2)e 2 + <3(e 3 ), with v 
the volume of a big sphere and e the small-to-large di- 
ameter ratio. It means that in the diluted regime of the 
small particles {P ~ p 2 ) this work can be estimated as 
- f3~ 1 {3/2)n 2 /e. In the" case of PHC AV = v e + O{e 2 ), 
with vo the volume of a big cube and e the small-to-large 
edge-length ratio. Again in the diluted regime of the 
small cubes the work is ~ /3 — x ?72 / e 2 . 




FIG. 4. Increment in the total volume volume, AV, avail- 
able to the small particles when two large particles come to 
touch each other (shaded region). If vo denotes the volume 
of a large particle, Av is (a) i>oe 2 (3/2 + e) for HS, and (b) 
v e(l + e 2 ), for PHC. 

In other words, the depletion induced by PHC is much 
stronger than that induced by HS. In the infinite asym- 
metry limit, the binary mixture HS has been shown to 
reduce to—the fluid of adhesive HS, provided 772 is kept 
constant.^ According to our estimation, in order to have 
a similar limit for the binary mixture of PHC we must 
scale the packing fraction of the small cubes as 772 = e£, 
with £ a constant. 

We can now assume this scaling of ?y 2 and take the 
e — > limit in Eq. (f42|). A tedious but straightforward 
calculation leads to 

c c ff (r) = cphcM + c ad (r) , (45) 

Cad(r) = [ {6 S (r) + yS(r) + 6y 2 V(r)} , (46) 

where y = r)i/(l — 771), cpHc( r ) is the DCF of the one- 
component PHC fluid [Eq. (56) of I], and 

5 s (r) = A(x, y, z) + A(y, z, x) + A(z, x, y) , (47) 
A(u, v, w) = S(ai - \u\)L(v)L(w) , (48) 
L(u) = (or - \u\)Q(ai - \u\) , (49) 

Q(t) being the usual Heaviside step function. Equation 
( p7| ) represents a delta function at contact of two large 
cubes, multiplied by the contact surface. The functions 
S(r) and V(r) are the overlap surface and volume, re- 
spectively, which already appear in the definition of cpnc 



[Eq. (53) of I] . In the zero density limit c e g becomes 

c eff (r)~/(r)+^ s (r), (50) 

with /(r) the Mayer function of the large cubes; so, as 
in HS, in the infinitely asymmetric mixture depletion in- 
duces an adhesive potential (in this case, of strength £). 
We will henceforth refer to this effective fluid as the fluid 
of parallel adhesive hard cubes (PAHC). 



C. Free energy functional of the fluid of parallel 
adhesive hard cubes 

Let us now take the e — > limit in the functional ( |30| ) 
to obtain the Hclmholtz free-energy functional for the 
fluid of PAHC. In the limit we will find that T — > 00; 
however, this is not a problem as long as for every fixed 
e there is a well-defined functional giving rise to a phase 
behaviour which do have a finite limit when e — > 0. We 
will show that this is the case, and thus this functional 
will be the effective functional we are looking for. We 
will see that the infinite contribution is just a constant 
shift in the origin of free energies, absolutely irrelevant 
for the phase behaviour. 

Let us begin by recalling what FMT prescribes for 
the semi-grand potential (|30|). It is convenient to in- 
troduce two dimensionless densities, r)(r) = afpi(r), and 
£(r) = crfe 2 /9 2 (r)- [There is no possible confusion be- 
tween the function r/(r) and the total packing fraction, 
because when e — > the total packing fraction is simply 
the packing fraction of the large component, i.e. the av- 
erage of n(r).] In what follows we will fix the unit length 
of our system by choosing oi = 1. In terms of these 
functions the FMT form of the functional ( |30| ) is 

/3T = d + J dr$(r) 

+e- 2 J dr C(r) (ln[V 2 e- 2 £(r)] - 1 - /3 M2 ) , (51) 
(3F id = J dr7 7 (r)(ln[V 1 7 ? (r)] - l) . (52) 

Equation (|5|) is just the ideal Hclmholtz free energy of 
the effective fluid, and <&(r) is given by (||), where, in the 
current notation, n a = 7/ * uj^ + e~ 2 u>^ . But £(r) is a 
dependent variable which should be eliminated in terms 
of 7y(r) and /i 2 via Eq. (|3l|), which in our case reads 

ln£(r)=ln*- V|l* w W(r), (53) 

a 

where we have defined the renormalised fugacity z = 
e 2 exp(/3/72 X ) ; f-T = 1^2 — /3~ 1 lnV 2 . For £(r) to have 
a well-defined expansion in powers of e we are forced to 
assume that z = 0(1). 
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We are almost ready to carry out the e-expansion. It 
only remains to determine the contributions to this ex- 
pansion coming from convolutions with 10% ■ Let f(u) 
be an arbitrary function of a single variable u. Then 

f*®= [ U+e/2 /(*) dt = ef{u) + 0(e 3 ) , (54) 

Ju—e/2 



f*6% = ±[f(u + e/2) + f(u-e/2)} 
= f(u) + ^f"(u) + 0(e i ). 



(55) 



Accordingly if /(r) is an arbitrary function of r, from the 
definitions (0) and the expansions above it follows 



f*4 a) = / + (e 2 /8)V 2 / + 0(e 4 ), 
f*ioi a) = e«/ + 0(e Q+2 ), 



(56a) 
(56b) 



for a — 3 or any vector component of a — 2 and 1. 
Then, assuming £ = £ + £i£ + £ 2 £ 2 + 0(e 3 ), the weighted 
densities n a can be expanded as 



n = £ e~ 2 + tie' 1 + |n + 6 + ^V 2 £ J 

ni = Coue- 1 + {£iu + Hi} + 0(e) , 
n 2 = {£ u + n 2 } + £ lU e + 0(e 2 ) , 
n 3 = n 3 + e£o + e 2 £i + e'^ 2 + 0(e 4 ) , 



where u = (1, 1, 1) and n a = r\ * w}' 



(a) 



0(e), (57a) 

(57b) 
(57c) 
(57d) 



From the expansions (|56]), fl57|), we can obtain 

,(<*) 



^^7* W2 



ln(l - 713) + — 



1 - n 3 



1. 



VH P 



1 - n 2 • n 2 



8£i+H x , (27/2)£ 2 + 4£ n 2 



.8' V 1 - "3 7 ' 2 (l-n 3 ) 2 l-n 3 (l-n 3 ) 2 
On the other hand, ln£ = In ft + §;£ + j|j - ^| e 2 + 0(e 3 ), hence (||) implies 

ft = ^(l - n 3 ) , 

! (l-n 3 ), 

37 



0(e 3 



£1 = -zn 2 - 8z 2 (l - n 3 ) 



6 = - o V 2 n 3 - 2W1 + 12z 2 n 2 + — z 3 (l - H 3 ) - $ ad , 



where for convenience we have introduced the shorthand 



. ld ^ z |Vn 3 | 2 -4n 2 -n 2 
8 1 - n 3 



(58) 



(59a) 
(59b) 

(59c) 



(60) 



We have already explicitly eliminated £ in terms of rj and 2, with the help of the expansion in e. We can now 
proceed to expand T itself, but before, let us rewrite the excess part appearing in ( |5l|) as 



In 



z(l - n 3 ) 



- 1 



+ f {$ + £- 2 £ln(l-n 3 )}; 



(61) 



then, using (|57j) and (|59j), and defining $ as in (g) with the n a replaced by n a (i.e., the excess free-energy functional 
of the one-component PHC fluid), 



+ e- 2 £m(l-n 3 ) = |zn 2 + 4z 2 (l - n 3 )} £ - 1 

119 



vr2 



z(l - n 3 ) 



9 -L -L £7 q . . Z Tit) Z 



£ = -z(l -n 3 ) + \ 8z 2 n 2 + 32z 3 (l - H 3 ) + - 



V • 



21 -H 3 



ln(l - n 3 )Vn 3 



0(e 3 



(62) 

$ + $ ad j + 0(e), 
(63) 







Substituting these two expansions into (|6l|), and using 
J V • [ln(l — n 3 )X7n 3 ] = 0, which holds if the density is 
constant at the boundaries or if it is a periodic function, 
we finally obtain 



T 



-U (e)V + ^(e)N + F + F ad + 0(e) , 



(64) 



where F is the FMT free-energy functional of the fluid 
of the large PHC; V and N are, respectively, the system 
volume and the number of large cubes; F ad — J $ ad is 
the new adhesive term; and Hq and /iq are given by 



55 



/?n (e) = ze~ 2 - 4zV x + — z 3 + 0(e) 



/3/Lt (e)=2e 2 + {3z - Az 2 }e~ 
+ t ^ _ I8z 2 + I 



0(e). 



(65) 



(66) 



The term — H n (e)V + po(e)N is divergent with e — > 0. 
It is the contribution of the small cubes to the free energy 
(as a matter of fact, their density is infinite in this limit). 
However it is irrelevant for the phase behaviour of the 
effective fluid because it simply adds Ilo to the pressure 
and fiQ to the chemical potential; as these two terms are 
independent of the density, they just cancel out in the 
equilibrium equations. Accordingly the final free energy 
functional for the effective one-component PAHC fluid 
turns out to be 

F PAHC ([phz)=F PW [p}+F ad ([p};z). (67) 

As a selfc onsi stency test, it is straightforward to show [us- 
ing Eq. (p9a|)] that c ad (r - r') = -f3S 2 F ad /5p(r)6p(r'), 



for the c ad function defined in fl46| ) 



D. Phase behaviour of the infinitely asymmetric 
binary mixture 

The phase behaviour of a very asymmetric binary mix- 
ture of PHC can be understood from that of the effec- 
tive fluid of PAHC, whose FMT free-energy functional 
we have just derived. 

As concerns the phase behaviour of the uniform PAHC 
fluid, from (|67|), (||), and ( |60| ) we can readily obtain the 
free energy per unit volume, 



/?/ = 77 {lnVx - 1 + In y + 3(1 - z/2)y + y 2 } 



(68) 



The pressure, P — —dF/dV = y 2 d(f /rj) /dy, turns out 
to be 



f3P = y + 3(l-z/2)y 2 + 2y 3 



(69) 



This equation has a van der Waals loop; the critical point 
can be found as the-,solution of the equations dP/drj = 
and d 2 P/drf = oJS i.e. 



1 + 6(1 - z c /2)y c + 6y 2 =0, 
1 - z c /2 + 2y c = 0, 



(70a) 
(70b) 



which is z c = 2(1 + y/2/3) ~ 3.63, and y c = l/v6, i.e. 

r] c = 1/(1 + V6) ~ 0.29. On the other hand, the equa- 
tion of the spinodal (dP/drj = 0) of this vapor-liquid 
transition is 

1 + 4t? + ?7 2 _ 
377(1-77) 

it is plotted in Fig. ||. Notice that this spinodal could 
have been obtained directly from (^) by taking the limit 
r — > 00, x — > 1, under the con strai nt 7*77(1 — x) = £ — > 
.2(1 — 77) [the limit follows from (59a)]. Thus it is not sur- 
prising that the line of critical points in Fig. || reaches 77 c 
for x — > 1. This makes clear the double interpretation of 
this transition: as a vapor-liquid transition of the PAHC 
fluid (with z^ 1 playing the role of a temperature), or as 
a demixing transition of the infinitely asymmetric binary 
mixture. 

Freezing of this system into a simple cubic lattice is 
again a continuous transition. Hence the transition line 
can be determine by the procedure described in Se c. , 
i.e. as the divergence of the structure factor [Eq. (jlT})] 
with the effective DCF (|45|) . The result is the line shown 
in Fig. ^. As it occurred for the general binary mixture 
(see Sec. |n|), the freezing line crosses the demixing spin- 
odal at a packing fraction smaller than 77 c ; in other words 
fluid-fluid demixing is a metastable transition. 

So far we have gone no further than we did in Sees. 
[II, IV. However this time we can study fluid-solid coex- 
istence because the density profile of the solvent is absent 
from the description. To proceed we again parametrise 
the density of the large cubes as in (||) , also with a pref- 
actor ■& to account for vacancies. We recall that the 
lattice parameter is related to this occupancy ratio by 
d = ('d/rf) 1 / 3 . The ideal contribution to the free energy 
per particle, \& , is again given by (14) (adding ln$ from 
the vacancies), and the hard-core part of the excess con- 
tribution is given by ( |l5| ) [with 713 (r) = r dp(x)p(y)p(z)\. 
We now need to work out the adhesive term. To this 
purpose first notice that 

00 

p'(u) 2 - q(uf = -4 g(u-nd+ 1/2) 

n.m— — 00 

xg(u- md- 1/2) ; (72) 
hence the adhesive free energy per particle can be written 
3 _22_ r d / 2 



dx 



d/2 



n,m= — 00 

xg(x- nd+ l/2)g(x- md - l/2)U(x) 

3 f°° 
— — ■dz I dx g(x) 



Y g(x-nd+l) 



n— — oo 



U(x + 1/2) , (73) 
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where U(x) is defined as 
U(x) = 



U(x) 



d/2 pd/2 

dy I dz 



p(y) 2 p{z) 



d/2 pd/2 

dyp'(y) / dzp'(z)V(x,y,z), (75) 

-d/2 J -d/2 



d/2 J -d/2 1 - $P(x)p(y)p(z) ' ( 74 ) V(r) = T(y)p(y)T(z)p(z) 



4-3n 3 (r)+3n 3 (r) 2 
1 - n 3 (r) 



(76) 



and it is periodic with period d [we have made use of 
periodicity in obtaining Q]. 



where T(u) = u if u G (—d/2, d/2) [the only relevant in- 
It is convenient to rewrite Eq. (|74|) integrating by part terval in J75|)] and it is d-periodic. Function V(x, y, z) is 



with respect to both variables, y and z; in doing so this 
equation becomes 



then also rf-periodic in all its three arguments; accord- 
ingly Eq. (fq) can be rewritten as 



U(x) 



and therefore 



ad 



dy g(y) / dz g(z) [V(x, y + 1/2, z + 1/2) - 2V(x, y + 1/2, z - 1/2) + V(x, y - 1/2, z - 1/2)] 



^2 g(x-nd+l) 



(77) 



■dz I dxg(x) / dyg(y) / dz g(z) 



n— — OQ 



x [V(x + 1/2, y + 1/2, z + 1/2) - 2V(x + 1/2, y + 1/2, z - 1/2) + V(x + 1/2, y - 1/2, z - 1/2)] , (78) 



also suitable for Gauss-Hermite numerical integration. 

In order to understand the effect of the adhesive con- 
tribution ([73|) let us see its asymptotic behaviour when 
a — > 00 and d — > 1 + (equivalently 1? — > ry + ), a limit which 
would represent a close packed solid. From its definition 
@ p{u) - 1 in this limit; thus U(x) - 1/(1 - rj). On 
the other hand, g(u) is very sharply peaked, so 

/oo 00 / q/ \ 1/2 

dxg(x) V fl (a; - nd + 1) ~ S (0) = (-) , 

therefore 

^~-|,-^p) 1/2 ^-oc. (79) 

2 1 — 77 V 7T / 

On the other hand 4" cx = 0(1) in this limit, while 
* id ~ 31ng(0) - (3/2) lna. In other words, the to- 
tal free energy per particle of the effective fluid mono- 
tonically decreases as the system approaches the close 
packing, regardless the value of density and solvent fu- 
gacity. This means that the system always collapses, i.e. 
the equilibrium phase behaviour is always a close-packed 
solid coexisting with an infinitely diluted gas. This sin- 
gular phase diagram is not exclusive of PAHC. For ad- 
hesive HS, the adhesiveness vs. packing fraction phase 
diagram [z plays the role of adhesiveness for PAHC) has 
recently been mapped out from simulations of the-square- 
well fluid in the limit of narrow and deep wells.cil These 
simulations prove that the only stable phases of this sys- 
tem are also a close-packed solid and an infinitely diluted 
gas. The reason fot,this pathology was put forward some 
years ago by Stelljia who showed that the partition func- 
tion of adhesive HS diverges if the number of particles 



is N > 12 (precisely the coordination number of the fee 
solid lattice). 





U 


MF 






MS \ 
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FIG. 5. Solvent fugacity, z, vs. solute packing fraction, 
r;, of the infinitely asymmetric binary mixture of PHC. The 
thick solid line separates the unstable region (U) from the 
metastable one; the thin one marks the (continuous) transi- 
tion from a metastable fluid (MF) to a metastable solid (MS); 
the dashed one is the fluid-fluid spinodal. 

In spite of the above we have seen that, as a function 
of a and the free energy per particle exhibits local 
minima at any value of z for some range of densities, 
the smaller z the wider this range. These local minima 
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correspond to metastable phases. The upper bound to 
the packing fractions at which local minima exist for a 
given z can be determined as the point where the com- 
pressibility vanishes. This upper bound, as a function of 
z, appears in Fig. |[ This figure shows the metastable 
phases. Notice that the region of metastability widens as 
z decreases. At low z these metastable phases are sepa- 
rated from the "collapse" by a large free energy barrier, 
so the system spends a long time in them before eventu- 
ally becoming a close-packed solid. As a matter of fact, 
if the system is prepared as a metastable solid at low z, 
for a long time it will show a pseudo-coexistence between 
two solid phases (an expanded solid and a close-packed 
solid). The situation resembles the isostructural solid- 
solid transition reported to occur in sorafecolloidal fluids 
with a narrow and deep attractive well.oH At higher z 
the same pseudo-coexistence should be observed between 
a diluted fluid and a close-packed solid. 

It is interesting to compare this phase behaviour with 
what has been determined to occur for adhesive HS using 
an effective-liquid DFT.ej Fluid and solid also appear as 
local minima of the free energy per particle as a function 
of the gaussian width; however, above a certain line the 
free energy becomes concave up to close packing. This 
puzzling behaviour was interpreted in Ref. ^5| as a per- 
colation transition. In the light of our findings it is the 
equivalent to the instability line of Fig. |[ The fluid-of ad- 
hesive HS also collapses into a close-packed solid.O The 
reason why this collapse has not been observed in Ref. 
[45| is that the theory used there does not account for va- 
cancies, and this forces the lattice parameter to be larger 
than 1 at any packing fraction. The instability manifests 
itself as the reported loss of convexity of the free energy. 



E. Polydispersity in the large cubes 

The singularity of the adhesive potential can be 
avoided b« .introducing polydispersity in the size of the 
particles E3'c£l It is clear that this prevents the system to 
form a perfectly packed solid. To see this effect on the bi- 
nary mixture we have introduced a small amount of poly- 
dispersity in the size of the large cubes. It is very easy 
to realise that starting off from a mixture of polydisperse 
large cubes and smal l cubes and repeating the process 
described in Sec. V C we end up with exactly the same 
form of the functional (j37|), with the n a 's now replaced 
by those corresponding to the polydisperse mixture. 

In order to make the simplest choice we consider the 
cubes as parallelepipeds and choose the length of each 
axis independently from a gaussian distribution of mean 
1 and variance Act. This particular choice has two im- 
portant advantages (they will be made clear below): (i) 
the free energy of the fluid phase is the same as that of 
the monodisperse system (hence its phase behaviour as 
well), and (ii) formally the expressions for the free en- 
ergy of the solid phase change very little. It also has two 



drawbacks: (i) particles are not cubic anymore, and (ii) 
there is a nonzero contribution in the negative lengths. 
As these two inconvenients disappear when Act — > they 
can be overcome by choosing Act <C 1. This choice also 
allows us to make two more simplifying assumptions: (i) 
the ordered phase must be a substitutional solid, i.e. the 
density profile can be expressed as p(r)V(cr), with V the 
normalised size distribution, er = (ct x , CTj,,.o,), and (ii) 
phase separation induced by polydispersityell can be ig- 
nored. 

Then, according to the definition of the n a 's 



n a (r) = J d(rV{a-)p*u)^ ) (r)=p*^ a ' ) (r), 



(80) 



i.e. it has the same definition as in the monodis- 
perse case, but the weights are redefined as uj^(r) = 
J derV(cr)uj^"\r). This amounts to replacing 6 U and 5 U 
in (0) by 



1 - erf ( V2 



>|-l/2 



Act 



S u = 



2ttAct 



{±} 



cxp < — 2 



>±l/2) 2 



Act 2 



(81a) 
(81b) 



which are like smoothed counterparts of the original 

weights. Since f°° du 9 u = f°° du5 u = 1. it follows 

j —00 j —00 

that the free energy of the uniform fluid is the same as 
that of the monodisperse system. Hence the fluid-fluid 
spinodal is the same as that shown in Fig. ||. 

We can determine the coexistence between the two 
fluid phases ,-by means of the usual double tangent 
construction.^! Figure ^ shows the resulting coexistence 
line. 

We can also assume a solid-like density profile as in 

the monodisperse case. Surprisingly enough, in spite of 

the striking difference of the smoothed weights defined 

by (|8l]) with respect to the original ones, when we obtain 

the corresponding weighted densities and work out the 

expressions a little bit, it turns out that the free-energy 

per particle of the polydisperse solid is simply given by 
_ ^poiy + yid + ^cx + ^ad^ w h ere the i ast t h ree con _ 

tributions are given by Eqs. (|l4|), (Jl5|), and (|7§|), with 
the slight modification that the parameter a appearing 
in the definitions ([H]) and (0) must be replaced by 



Q 



1 + qAct 2 /2 



(82) 



(of course, in these expressions a = 1, the mean value, 
and \& ld carries the additional In 1) to account for vacan- 
cies), and where ^ poly = — In (\Z2~7tAct) — 1 is the entropy 
of mixing (an irrelevant constant). 

From Eq. (|8^) it can be seen that no matter how small 
Act be, for small a's (a <C Act -2 ) a ~ a, and the sys- 
tem is "blind" to polydispersity, whereas for large a's 
(a » Act -2 ) a ~ 2Act -2 , i.e. the system never collapses. 
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As a consequence the singular behaviour of the monodis- 
perse system is removed, and we can readily determine 
phase equilibria. A typical result for a small value of Act 
is shown in Fig. |[ This figure reveals several remarkable 
features. Firstly, it shows that the fluid-fluid transition 
is metastable. Secondly, there is an isostructural solid- 
solid transition for a certain interval of z. In this inter- 
val the expanded solid (Si) appears after a continuous 
transition from the fluid phase (F). The expanded solid 
(Si)-dense solid (S2) transition ends at a critical point, 
z s , below which we can only find a fluid and a single 
solid, separated by a continuous transition. Thirdly, as 
z increases from z s the expanded solid packing fraction 
decreases down to meeting the freezing packing fraction. 
Above the point where this occurs coexistence is between 
a fluid and a dense solid (S2), the former quickly be- 
coming highly diluted and the latter highly packed. No- 
tice the strong resemblance between this true equilibrium 
phase behaviour and the metastable behavio ur de scribed 
to occur the monodisperse PAHC fluid (Sec. |V D| ) . 



N 




0.5 1 

FIG. 6. Same as Fig. ^ for the polydisperse PAHC fluid 
(Act = 0.045). The thick solid line marks the fluid-solid or 
solid-solid coexistence; the thin one marks again the (contin- 
uous) fluid-solid transition below the coexistence region; the 
dotted line is the metastable fluid-fluid coexistence. 



difficult (the closely related fluid of parallel hard paral- 
lelepipeds has been recently used in— the context of as- 
sociating fluidscS and liquid crystalsea). The main con- 
tribution of the fluid of PHC is to the understanding 
of the phase diagram of a binary mixture. TJbif fluid 
proved to undergo stronger depletion than HSBEl How- 
ever, as it has been shown in this work, this feature is 
irrelevant when the effect of depletion in the solid phase 
is accounted for. Spatial order of the large component 
strongly enhances demixing, so that fluid-solid demixing 
becomes the main scenario of the phase diagram of bi- 
nary mixtures. But this transition can be preempted by 
the freezing of the large component, and when this hap- 
pens the system phase separates into two fee solids with 
a different lattice parameter. This effect, very clearly 
shown here for the mixture of PHC (in the limit of infi- 
nite asymmetry), has also been confirmed in simula tions , 
of HS interacting via an effective depletion potentialj23S 
and very recently also in direct simulations on the true bi- 
nary mixture.E3 The simulations also show that the solid- 
solid transition disappears as the asymmetry of the two 
components decreases, but it is anyhow categoric with 
respect to the fluid-solid nature of demixing. 

A final remark concerns the two-dimensional mixture. 
We have made preliminary calculations in this case and 
have found an adhesive contribution similar to the three- 
dimensional one. We have nor carried out a detailed anal- 
ysis yet, but the same collapse is present in this case, thus 
indicating a behaviour qualitatively similar to the one 
shown here, except that we cannot say anything on the 
existence of a solid-solid transition. This results arejia 
perfect qualitative agreement with recent simulations.Ej 
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VI. DISCUSSION AND CONCLUSIONS 

The fluid of PHC is a rather academic one which 
however has the great advantage of being analytically 
tractable in contexts where the fluid of HS is not, thanks 
to its adequacy to a fundamental measure description. 
Yet, with some .peculiarities due to the lack of rota- 
tional symmetryElu the physics it reveals is similar to 
that of more realistic fluids. It then allows for theoreti- 
cal investigation on fluid phase behaviour otherwise very 
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